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Abstract 

In the context of the finite elasticity theory we consider a model for 
compressible solids called "compressible neo-Hookean material". We 
show how finite-amplitude inhomogeneous plane wave solutions and 
finite-amplitude unattenuated solutions can combine to form a finite- 
amplitude Love wave. We take a layer of finite thickness overlying a 
solid half-space, both made of different pre-stressed compressible neo- 
Hookean materials. We derive an exact solution of the equations of 
motion and boundary conditions, and also obtain results for the energy 
density and the energy flux of the waves. Finally, we investigate the 
special case when the interface between the layer and the substrate is 
in a principal plane of the pre-stain. A numerical example is given. 

1 Introduction 

A seismic event launches at least two types of surface waves, one causing ver- 
tical (elliptic) movements, the other causing rather destructive lateral (shear 
horizontal) movements. In the essay which won him the Adams Prize in 1911, 
Love [1] proposed a simple Earth model which supports the latter kind of 
waves, by considering a crust made of an isotropic, linear, elastic solid, rigidly 
bonded onto a substrate (a semi-infinite solid) made of another isotropic, lin- 
ear, elastic solid. In this heterogeneous structure, a shear horizontal wave 
may propagate, leaving the upper face of the layer free of traction, and hav- 
ing an amplitude which decays rapidly with depth in the substrate. This 
localisation of the amplitude variation is what makes Love waves (and sur- 
face waves in general) a subject of great interest in seismology, because it 
means that the energy spreads essentially in two dimensions, and thus the 
wave travels further from the epicentre than bulk waves, for which the energy 
spreads in three dimensions. 
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Over the years, Love's results were extended in several directions and 
turned out to be also useful in other contexts. For instance anisotropy (e.g. 
[2]), inhomogeneity (e.g. [3], HI Ej), piezoelectric coupling (e.g. [SI [TJ [8]), 
and many other effects were considered — see the review by Maugin [9] for 
an exhaustive account. Two possibilities are of special interest to seismol- 
ogy science: the possibility of including strain-induced anisotropy, because 
it provides a simple and revealing modelling of the consequence of slow tec- 
tonic movements, and the possibility of including non-linear effects, because 
large amplitude seismic movements have indeed been observed. The first 
possibility can be dealt with within the framework of small-amplitude waves 
superimposed upon a large static homogeneous pre-strain, see Hayes and 
Rivlin [10], Willson [H], Kar and Pal [12], or Dowaikh [IS]. The second 
possibility is usually treated in the framework of the so-called weakly non- 
linear elasticity theory (see Maugin [11] or Norris [15] for instance), where the 
equations of motion are developed one order further than the linear regime 
(Note that according to Zabolotskaya [16], non-linear shear horizontal waves 
require fourth-order elasticity.) 

Here we propose to combine both effects, by considering finite-amplitude 
Love waves in a finitely deformed layer/substrate structure. Our only re- 
striction is in our choice of a constitutive law for the solids, as we focus on 
the so-called compressible neo-Hookean materials. These enjoy a peculiar 
property once deformed [T71 [THl [IH] : they allow the propagation of a finite- 
amplitude, inhomogeneous, linearly polarized, transverse plane wave in any 
direction. Moreover, this wave is obtained by solving a linear ordinary differ- 
ential equation even though the theory is completely non-linear (In passing 
we note that the eventuality of a solitary wave is thereby precluded here). We 
use this wave as an ingredient in the construction of the Love wave solution 
to the corresponding boundary value problem. A great deal of generality is 
nonetheless achieved, in particular because non-principal wave propagation 
is possible, and because the solution is exact, without any limitations to be 
imposed on its magnitude. Corresponding general results are obtained on 
energy propagation in the layer and the substrate. 

A special case of the 'compressible neo-Hookean' model was first intro- 
duced to describe a class of solid polyurethane rubbers studied in the Blatz- 
Ko experiments [28]. Of course we remain aware that solids playing a role 
in the applications of Love waves are not always adequately described by 
the 'compressible neo-Hookean' strain energy density. However exact and 
expilicit results like those of this paper are always useful, either as a basis 
for perturbation methods or for testing numerical schemes. 

The paper is organised as follows. In Section [2], we present the consti- 
tutive equations of the compressible neo-Hookean model. Next (Section [3]), 
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we retrieve results for transverse inhomogeneous time-harmonic waves su- 
perimposed on a pre-stressed state [TTl [TS] but also present similar results 
for transverse unattenuated time-harmonic motion. In Section HI the set-up 
consisting of a semi-infinite substrate covered with a layer of finite thickness 
is described and the effect of pre-strain is investigated, assuming that the two 
solids are rigidly bonded. A large amplitude Love wave solution is then ob- 
tained provided the propagation direction in the interface and the normal to 
the interface are along conjugate directions of the pre-strain tensors. A dis- 
persion equation is also obtained, similar to that of linear isotropic elasticity. 
Here however this equation involves the bulk wave speeds along the propa- 
gation direction and along the normal to the layer, which are both affected 
by the pre-deformations. In Section [Ul we derive the energy densities and 
energy fluxes corresponding to the motions in the layer and in the substrate. 
Mean energy densities and fluxes are obtained by averaging over a period in 
time, and their properties are investigated. Total energy densities and fluxes 
are also introduced, showing that the repartition of energy between the layer 
and the substrate depends on the ratio of layer thickness to the wavelength. 
Finally (Section [7j), we consider the special case when the interface is in a 
principal plane of the pre-strain tensors in the layer and the substrate. In 
this case the propagation direction may be any direction in the interface and, 
owing to the anisotropy induced by the pre-stress, the Love wave speed varies 
with the propagation direction. 



2 Compressible neo-Hookean materials 

In order to deal with possibly large deformations of solids, we invoke the 
finite elasticity theory. To model the non-linear elasticity of solids, we use 
one of the simplest constitutive models on offer for compressible isotropic 
solids, which may be called the 'compressible neo-Hookean material'. We 
here present this model. 

Let F be the deformation gradient, defined as usual (see, for instance, 
EH) by 

F = dx/dX, FiA = dxi/dXA, (2.1) 

where X is the position of a particle in the reference (Lagrangian) configura- 
tion and X the corresponding position in the current (Eulerian) configuration. 
Associated with F is the the left Cauchy-Green strain tensor 

B = FF^, Bij = {dx^/dXA)idxj/dXA), (2.2) 

where the superscript T denotes the transpose. The determinant of F is 

J = det F = {det By/\ (2.3) 
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It is the ratio between the volume of a material element of the solid in the 
reference and current states. Then, compressible neo-Hookean materials are 
characterized by a strain-energy density W, measured per unit volume in the 
undeformed state, given by 

2W = fi{tTB -3) + G(J) -G(l), (2.4) 

where is a constant (the shear modulus) and G{J) is an arbitrary function of 
J (a material function, which can be adjusted to model the compressibility 
properties of the solid). The corresponding constitutive equation for the 
symmetric Cauchy stress tensor T is 

T = lG'{J)I + fiJ^^B. (2.5) 

Such a 'compressible neo-Hookean' material is sometimes called 'special Blatz- 
Ko' material, or 'restricted Hadamard' material [T7 1 fT8l [T9] . Hayes [20] shows 
that the resulting equations of motion are strongly elliptic when 

/i > , G"{J) ^ 0, (2.6) 

and we assume as much henceforth. We also assume that the undeformed 
state is stress free, so that G'{1) = —2^. Note that comparison with linearised 
isotropic elasticity yields G"{1) = 2(A -|- fi), where A and /i are the Lame 
coefficients. Here however, no restriction is placed on the amplitude of the 
displacement u = x — X, so that due to (12.21) . and the arbitrariness of G{J), 
the relation between T and u is clearly non-linear. 

Several examples of specific volumetric functions G{J) have been pre- 
sented over the years (see, for instance, Ba§ar and Weichert [21]). Among 
these, we recall the particularly simple choice of Levinson and Burgess [22] . 
leading to a model called 'simplified Blatz-Ko' in [TH], that is 

G(J) = (A + /i)(j2 - 1) - 2(A + 2/i)(J - 1), (2.7) 

with /i > and A -|- yU > in order to satisfy (12. 6p . 

3 Transverse waves superimposed on a static 
deformation 

Here we consider transverse wave solutions in unbounded pre-strained ma- 
terials. Suppose that a compressible neo-Hookean material is first subjected 
to a static finite homogeneous deformation defined by 

X = FX, Xi = FiAXA, (3.1) 
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where the Fia are constants. The corresponding constant left Cauchy-Green 
strain tensor is B = FF"'" and the determinant J = det is a constant. 
On this state of deformation, we superpose a time-dependent displacement 
taking a particle from position x to position 

x = x{x,t) = x{FX,t) = x + u{x,t), (3.2) 

where u is the mechanical displacement. In the absence of body forces, the 
equations of motion for this time-dependent deformation may be written in 
the form [HI [IS] 

px = div^P, pxi = dPik/dxk, (3.3) 

where p = PqJ~^ is the constant mass density in the intermediate state of 
static deformation, and P is the Piola-Kirchhoff stress tensor at time t with 
respect to the intermediate state of static deformation, 

P = (det F)TF~^, = (det F)T,,F^/. (3.4) 

Here T is the Cauchy stress tensor at time t, and F = dx/dx is the defor- 
mation gradient with respect to the intermediate state of static deformation. 
We note that F = FF~'^, where F = dx/dX is the deformation gradient 
with respect to the undeformed configuration. 

As in [IH], we are now looking for solutions with a displacement field 
u = x — X of the form 

u = f{m ■ x)g{n ■ x — vt)a, (3.5) 

where / and g are functions to be determined, and m, n and a are unit 
vectors. It is assumed that m and n are not parallel and that a is orthogonal 
to both m and n, 

a ■ m = a ■ n = 0, (3.6) 

so that fl3.5p represents a linearly polarized transverse wave with propagation 
speed V. 

For such a wave motion, recall that F and its inverse are [T8] 

F = I + f'ga®m + fg'a®n, F ^ = I - f ga® m - f g' a® n, (3.7) 
so that the special Blatz-Ko constitutive equation (12. 5p yields 

T = \G\J)I + ^iJ-^FBF^ , (3.8) 
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and, because det = 1, the Piola-Kirchhoff stress tensor (13.41) reduces here 
to _ _ 

P = lG'{J)F ^ + nJ'^FB. (3.9) 

It follows [18j that the displacement field (13. 5p is a solution of the equa- 
tions of motion if and only if / and g satisfy the equation 

{n - Bn- fi'^poV^)fg" + 2n ■ Bmf'g' + m ■ Bmf'g = 0, (3.10) 

where /' and g' denote the derivatives of / and g with respect to their 
argument. We then choose m and n such that [17\ 

nBm = 0, (3.11) 

which means that m and n are along the principal axes fo the elliptical 
section of the ellipsoid x • Bx = 1 by the plane a • x = 0. Then, (I3.10p 
yields two uncoupled equations for / and g : 

vL{f"/f) = -c, {vl-v'){g"/g) = c, (3.12) 

where c is an arbitrary constant, and Vm and Vn are the wave speeds of 
homogeneous bulk waves propagating along m and n, respectively, 

Po'^m = prn- Bm, pov"^ = pn- Bn. (3.13) 

If c is assumed to be positive, c = k"^ (say) for some real k, then (I3.12p 
yields an unattenuated time-harmonic wave motion provided f ^ > f^. The 
displacement field of this wave is 



u{x, t) 



B sin( — m ■ x) + C cos( — m ■ x) 



V w - vi 

where k, B, C are arbitrary constants. 

If c is assumed to be negative, c = —7^ (say) for some real 7, then (I3.12p 
yields an inhomogeneous time-harmonic wave motion provided f ^ < f^. The 
displacement field of this inhomogeneous plane wave is 

u{x, t) = Aexp ( — —m ■ x ) cos ( — ^ {n. ■ x — vt) la, (3.15) 



where 7, A are arbitrary constants. Here we retrieve a solution obtained in 
p!7l lis] . However the solution (I3.14p was not mentionned in these papers 
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because the emphasis there was on inhomogeneous plane waves. Here, both 
(I3.14P and (13.151) are needed for the construction of a Love wave solution. 

We remark that when the condition (13. lip is not satisfied, solutions may 
nevertheless be obtained [18]; however either f ot g is then of real exponential 
type and hence no time-harmonic wave motion is possible. In this paper we 
focus on time-harmonic waves because they are the building blocks for Love 
waves. 

For future reference, we conclude this section with the evaluation of the 
traction vector t on a plane m ■ x = constant . Because the displacement u 
is along a and hence orthogonal to m, such a plane is globally preserved in 

the motion. Moreover, det = 1 and F m = m by (\3.7h . Hence, using 

Nanson's formula, da = (det F)F da, linking an areal element da in the 
current configuration with the same areal element da in the intermediate 
state, we conclude that when da is along m, then the areal element is the 
same in both configurations: da = da = mda. Thus, the traction vector t 
on a plane m ■ x = constant is the same whether it is measured per unit 
area of the intermediate state or per unit area of the current state. Recalling 
(13. lip , we obtain 

t = Pm = Tm = Tm + pv'i,fga, (3.16) 
where T denotes the constant Cauchy stress tensor of the intermediate state. 

4 The pre-stressed layered formation 

We wish to extend the classical results of Love [Ij in the linear elasticity the- 
ory in two directions: by taking account of initial stresses (and the accom- 
panying strain- induced anisotropy), and by allowing the wave's amplitude to 
be arbitrarily large. 

We start with Love's original set-up, which consists of a semi-infinite 
substrate, covered with a layer of finite thickness. The two solids are bonded 
rigidly. Here we assume that both the substrate and the layer are made of 
different 'compressible neo-Hookean materials', with a shear modulus /i and 
a function G for the substrate and a shear modulus Ji and a function G for 
the layer. Also, po and po denote the mass densities of the substrate and the 
layer, respectively, measured in the undeformed reference configuration. 

In order to model geological formations, it is common to consider that 
the solids have been subjected to initial stresses, giving rise to strain-induced 
anisotropy (see for instance the works of Biot [23j or Tolstoy |24j). To simplify 
matters here, we focus on static homogeneous initial strains. Thus, if X 
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denotes the position of a material particle in the undeformed solids, with 
origin X = in the interface between the substrate and the layer, then the 
initial deformations are 

X = FX, X = FX. (4.1) 

Here, the components of the deformation gradients F in the substrate, and 
F in the layer, are constants. The associated constant left Cauchy-Green 

strain tensors are B = FF in the substrate and B = FF in the layer. 
Also, we let J = det F and J = det F. 

We call m the unit vector normal to the faces of the layer, in the static 
pre-strained state (14.11) . oriented from the layer toward the substrate, see 
FigUl Hence, in this pre-strained state, the layer /substrate interface is the 
plane m • a; = 0, or equivalently, m ■ x = 0, and the substrate occupies 
the m • X ^ half-space. Also, in the static pre-strained state, the upper 
face of the layer, in contact with vacuum, is the plane m ■ x = —h, where 
h is the thickness of the layer in this state, so that the layer occupies the 
—h < m ■ X < region. 

As in the classical linear case, we focus on the possible existence of a 
linearly-polarized transverse wave, propagating in a direction n and polarized 
in a transverse direction a, both parallel to the interface. Thus, (n, a, m) 
forms an orthonormal triad. 

To ensure rigid bonding the displacement must be continuous at the in- 
terface. The displacements u in the substrate and u in the layer are given 
respectively by 

u = X - X = {I - F-^)x, u = x- X = {I -F~^)x. (4.2) 

At any point in the interface m • x = m ■ x = 0, we have x = x = an + (3a 
for some a and (3. Then the displacement continuity u = u a.t the interface 

is equivalent to F~^{an + Pa) = F (an + Pa). Because this must hold for 
all a, P, we conclude that the requirement of displacement continuity is 

F^n = F~^n, F^a = F^a. (4.3) 

— 1 ~-i 

It then follows that -F nxF a = F nxF a. Hence, using the identity 
[25]: F-^n x F'^a = {det F)-^F'^{ n X a), and similarly for F, we obtain 

(det F)-^F'^m = (det F)-^F^m, (4.4) 
from which it follows that 

J^^m ■ Bm = J^^m ■ Bm. (4.5) 
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m.x=-h 



m .X = 



m.x>0 




Figure 1: Sketch of the layered formation. The wave propagates in the 
direction of the unit vector n, it is polarized along a, and the magnitude of 
its amplitude varies along the direction of the unit vector m, normal to the 
faces of the layer. 

We note in passing that the unit vector M normal to the faces of the layer 
in the undeformed state is given by 

M = (m ■ Bm)-^/'^F'^m = (m ■ BmY^''^F^m, (4.6) 

and that the thickness H of the layer in the undeformed state is 

H = h{m- Bm)-^/'^. (4.7) 

Next, we consider the corresponding constant Cauchy stress tensors T 
and T, in the substrate and the layer, respectively. The equihbrium of the 
pre-strained state requires that the upper face of the layer be subjected to 
the traction (deadload) r = Tm, and that the traction vector be continuous 
at the interface, Tm = Tm. Using the constitutive equations of the layer 
and of the substrate, this yields 

JlJ^^n ■ Bm = fiJ^^n • Bm, JlJ^^a ■ Bm = jiJ^^a • Bm, 
\G'{J) + JiJ-^m ■ Bm = \G'{J) + ^J-^m ■ Bm. (4.8) 

Using the requirements of continuity of the displacement and of the trac- 
tion at the interface, we now show how a given initial strain B in the layer 
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(resulting from a prescribed stress T) determines the initial strain B in the 
substrate. 

First, we note that (jl!8])3 may alternatively be written as follows, using 

iG"( J) + fiJJ-^m ■ Bm = iG"( J) + JlJ-^m ■ Bm. (4.9) 

Because B is given, J and m ■ Bm are known and so, this is an equa- 
tion for a single unknown, J. It is shown in the Appendix how the strong 
ellipticity assumption implies that this equation has at most one positive 
solution for J, and one and only one solution if, in addition, it is assumed 
that limj_»o G'{J) = — oo. Once J is known, then (14. 5 p determines m • Bm, 
and (14.81) 19 determine n • Bm and a • Bm. Explicitly, 

m • Bm = (J/jYm ■ Bm, n • Bm = (J/ J)(jl/ fi)n ■ Bm, 

a - Bm= {J/J){Jl/i2)a-Bm. (4.10) 

In order to obtain all the components of B in the orthonormal triad 
{n,a,m), we still need to determine a • Ba, n • Bn, n • Ba. For this 
purpose, we note that the displacement continuity requirement (14.31) implies 
that 



n ■ B n = n ■ B n, a ■ B a = a ■ B a, n ■ B a = n ■ B a. 

(4.11) 

Using n = a X m in (I4.1ip i . and the identity 

n ■ B = J"^[(a ■ Ba){m ■ Bm) - {a ■ Bmf], (4.12) 

valid also with B instead of B, we obtain 

„ ~ (a ■ BmY {a ■ Bm)^ 

a Ba = a Ba~^ + ^ -. 4.13 

m ■ Bm Bm 

Recalling (l4.10P i we obtain 

^ fji^ \ (a- Bmf 

a- Ba = a- Ba+ [^-l]^ (4.14) 

J m ■ Bm 

A similar procedure, using a = m x n in (14.110 9. yields 

„ ^r; fT^^ \ in ■ BmY 
n- Bn = n- Bn+ -l]^ (4.15) 

V/" / m ■ Bm 
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Finally, (14.1 II) ';; with n = a x m and a = m x n yields 



nBa = nBa+{^-l 

/i2 



(n ■ Bm){a ■ Bm) 
m ■ Bm 



(4.16) 



To summarize: when the (constant) left Cauchy-Green strain tensor B in 
the layer is prescribed, then the (constant) left Cauchy-Green strain tensor 
B in the substrate is uniquely determined. First, J is uniquely determined 
from equation (14. 9p . Then, all the components of B in the orthonormal triad 
{n,m,a) are explicitly given by equations (I4.10l) i 9 •;!■ 04.141) . fl4.15p . (14.161) . 

In order to use the exact wave solutions described in Section [HI we shall 
assume from now on that the condition (13. lip is fulfilled in the layer, which, 
by (14.1 op 9 implies that it is also fulfilled in the substrate : 



n ■ Bm = 0, n ■ Bm = 0. 

Then, the equations for the other components of B reduce to 



(4.17) 



m • Bm = (J/J^m ■ Bm, a • Bm = {J/ J)(n/ fi)a ■ Bm, 

■2 ^ (a BmV 



a ■ Ba = a Ba + ( ^ - 1 



m ■ Bm 



n ■ Bn = n ■ Bn, 



n ■ Ba = n ■ Ba. 

In tensorial form, the expression of B in terms of B is given by 
(m ■ Bm)B = (m ■ Bm)B — Bm Bm 



(4.18) 



+ 



{jj,/ li){a ■ Bm)a + {J / J){m ■ Bm)m 

{Jl/lj){a-Bm)a+{J/J){m-Bm)m . (4.19) 



In particular, we note that Bn = Bn. 



5 Large amplitude Love wave 

Now we look at wave propagation in the initially deformed structure. In the 
substrate, we require that amplitude of the wave decays in the direction of 
m, and hence the displacement field u is assumed to be of the form (13.150 . 
In the layer, we consider an unattenuated time-harmonic displacement field 
u of the form (I3.14p . As in the classical case, we wish to combine these exact 
wave solutions in order to obtain a global time-harmonic wave motion with 
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propagation speed v. Note that both displacement fields need to be of the 
same angular frequency, and hence the same wavenumber, in order to satisfy 
boundary conditions at the interface. Thus, using (13.151) for the substrate 
and (13.141) for the layer, we write 

u{x,t) = Aexp ( — —m ■ x ) cos k{n ■ x — vt)a, (5.1) 



m 



and 

u{x, t) 
where 



i?sin(::^m ■ x) + C cos{:^m ■ x) 



K 7 



V 



COS k{n ■ X — vt)a, (5.2) 



(5.3) 



2 ' 



is the wavenumber. Here, in accordance with (13.131) . the body waves speeds 
Vn, Vm, Vn, Vm are giveu by 

pov'^ = pm • Bm, p^v"^ = pn • Bn, 



n 

PoV^ = P^- Bm, povl = Jin- Bn, (5.4) 

and the Love wave speed v has to satisfy 

< < (5.5) 

Notice that this is possible only when v"^ > v^, or equivalently, recalling 
n • Bn = n ■ Bn, when fi/po > p/po- Thus, Love waves require the combi- 
nation of a 'slow' (or 'soft') layer over a 'fast' (or 'hard') substrate, indepen- 
dently of the initial pre-strain. 

We now show that the boundary conditions may be satisfied. This leads 
to the dispersion equation (a relation between k and v) and the determination 
of the constants A, B, C in terms of a single parameter characterizing the 
amplitude of the wave. 

The first boundary condition to enforce is that the displacement is con- 
tinuous at the layer /substrate interface m • x = 0. Using (15. ip and (15.21) . 
this gives 

A = C. (5.6) 

The second boundary condition is the continuity of the traction vector at 
the layer /substrate interface m ■ x = 0. Using (I3.16p . and applying it to the 
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wave motion flS.ip . we obtain, for the traction t on a plane m - x = constant 
in the substrate, 



t = Tm — 7pf^y4exp j — —m ■ x j cos k{n ■ x — vt)a, (5.7) 

where p = J~^po is the mass density of the substrate in the intermediate 
configuration. Similarly for the traction i on a plane m ■ x = constant in 
the layer, we find 



t = Tm+Kpvj 



B cos I ^:^m ■ X 

Vrr,. 



C sin I ■ x 



COS k{n ■ X — vt) a. 

(5.8) 

where p = J ^po is the mass density of the layer in the intermediate config- 
uration. Also, recall that Tm = Tm = r, where r is the constant traction 
(deadload) applied at the upper face of the layer. Hence the condition t = t 
at the layer /substrate interface m ■ x = reads 



jpVmA + upVmB = 0. 



(5.9) 



The third boundary condition is that, at the upper face of the layer 
m ■ X = —h, the wave creates no traction in addition to the static traction 
(deadload) r. Using fl5.8l) . this yields 



B cos{:^h) + C sm{:^h) = 0. 



(5.10) 



Equations (15. 6p . (15.91) . (15.101) form an algebraic linear homogeneous sys- 
tem for the three unknowns A, B, C. Writing the condition for non-trivial 
solutions and using (15.31) . we arrive at the following dispersion equation re- 
lating the wave speed v to the wave number k, 



tan 



kh\ 



I y2 _ y2 



pVr 



pVr 



(5.11) 



Let c and c denote the transverse bulk wave speeds in the underformed 
substrate and layer, respectively: = /u/po, = p/'po- Using (I4.18P and 
(15. 4p . we note that v^/v^ = /T? and pvm/ipvm) = Poc/iftoc), so that the 
dispersion equation (15. lip may also be written as 



tan 



PqC [cjc)^ - {v/Vnf 
PoCy {v/VnY - 1 



(5.12) 
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20 25 

(v^vjkh 



Figure 2: Dispersion curves for successive Love wave modes in a pre-strained 
configuration. Here Hi = /7r/y^(c7cp"^-T, (/ = 1,2,3,4,...). Thus, for 
< kh < (vm/vn)Hi, only one mode may propagate (fundamental mode). 
For {vm/vn)Hi < kh < (ym/vn)Hi^i, I + 1 modes may propagate. 



In the absence of pre-strain, B 
in the substrate and 
equation specializes to 



IP 

m 



- B = I, J = J = 1, hence f ^ = f ^ = 
in the layer, so that this dispersion 



tan 



kh^y (f /c)^ 




c/c)2 - (f/c)^ 
(t;/B)2 - 1 



(5.13) 



which coincides with the dispersion equation of linear isotropic elasticity [29] . 
From a practical point of view, any dispersion curve obtained from fl5.13p as 
a plot of v/c against kh for a given choice of po/Po and c/c can be used 
in the present context of f l5.12p . by identifying f/c with v/vn and kh with 
{vn/vm)kh, see Willson [llj for similar results in the small-on-large theory. 
A typical plot for the different wave modes is presented in Figj2l Of course, 
the scope of the results is now richer because they include large amplitudes 
and pre-stress. 

When the dispersion equation (15. lip is satisfied, the solution of the linear 
homogeneous system (15. 6p . (15. 9p . (I5.10p for A, B, C is 



A = a cos ( ::^h 



B 



—a sin ( ■:^h 



C = a cos I ::—h 



(5.14) 
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where a is arbitrary. Thus A, B, C are expressed in terms of a single pa- 
rameter a characterizing the amphtude of the Love wave. 

Finally, we now denote by [r], ^, () the Cartesian coordinates along (n, a, m) 

1] = n ■ X, ^ = a ■ X, ( = m ■ X, (5.15) 

and we find that the displacement fields (15. ip (15.21) are 

u = a cos ( :z:—h ) exp ( — —( ] cos k(ri — vt)a, 

\Vrr,. J \ Vm. J 



K 



u = a cos ^:^{h + COS kir] — vt)a, 
or equivalently, recalling (15.31) . 



(5.16) 



u = a cos 



u = a cos 



khi 



exp 



k{h + o, 



' f ^ — f ^ 



COS k{ri — vt)a. 



cos k{ri — vt)a, (5.17) 
(5.18) 



6 Energy density and energy flux 

Here we compute the energy flux and the energy density associated with a 
motion of the type (13.21) in compressible neo-Hookean materials. Let W and 
W be the strain energy densities corresponding, respectively, to the static 
deformation (13. ip and to the motion (13. 2p . both measured per unit volume 
of the undeformed state. Then, the energy density S, measured per unit 
volume of the homogeneously deformed state (13. ip . and the corresponding 
energy flux vector TZ are [27J 

£ = lpW-W+J-\W-W), nu = -%'Pik, (6.1) 

where the Piola-KirchhofF stress tensor P with respect to the state of ho- 
mogeneous static deformation is defined by (13. 4p . They satisfy the energy 
balance equation 

{dS/dt) + {dTlk/dxk) = 0, (6.2) 

where Xk are the coordinates in the state of homogeneous static deformation, 
and the partial derivative with respect to time is taken at fixed x. 

We now evaluate the energy density and the energy flux vector for the 
wave motion (I5.I6P 9 in the layer, and for the wave motion (15. 160 1 in the 
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substrate. Using (15 .161) 9. we obtain for the layer 

I = lpk^a'^{v'^ + ^) cos2 + sin^ k{r] - vt) 

+ \'pk^a^{v^ - IPJ sin^ ^{h + C) cos^ k{7] - vt) (6.3) 
—a a- sin 7:^{h + Q cos k{ri — vt)Tm 

+k cos (/i + C) sin k{ri — vt)Tn], 

Vm 

72. = cos ^:^{h + Q sin k{r] — vt){—vkaTa (6.4) 
+a^kvJlJ^^[:z:— sin ^::—{h + Q cosk{rj — vt)Bm 

^ , , 

+k cos (/i + C) sin k{ri — vt)Bn]}. 

Using fl5.16P i . we obtain for the substrate 
8 = \pk'^a'^ cos'^{^h)exp{^() [v^ + v"^ sin^ k{r] — vt) — v"^ cos^ k{r] — vt)] 

—a cos(:c— /i)exp( — () a- [ — cos k{r] — vt)Tm 

+k sin k{r] — vt)Tn], (6.5) 
72. = —vka cos{^:^h) sin k{r] — vt)exp{ — —QTa (6.6) 

+a'^kvjj,J~^ cos'^{^:^h)exp( () [ — cos k{r] — vt)Bm 

+k sin k{r] — vt)Bn] . 

Mean energy densities and mean energy fluxes are obtained by averaging 
over a period in time at fixed x, 

{£) = {uj/2n) S{x,t)dt, {n) = {uj/2n) n{x,t)dt. (6.7) 

Jo Jo 

Using (16. 3p . (16. 4p . (16. 7p . we find the mean energy density and the mean 
energy flux in the layer as 

(8) = \pk^a\v^ + ?^[cos2 + - sin^ + 0]}, 

{Ti) = \vk'^a^JlJ-^ cos2 #-(/i + C) Bn. (6.8) 

Using (16. 5p . (16. 6p . (16. 7p . we find the mean energy density and the mean 
energy flux in the substrate as 

{S) = \pviea^cos\J^h)eM^O. 

{H) = \vk^a^pJ-^ cos^^h) exp(^C) Bn. (6.9) 
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The mean energy flux in the layer and the mean energy flux in the substrate 
are both along Bn = Bn, thus along the same direction, parallel to the 
interface. In general, this direction is not along the propagation direction n. 
This is an effect of the anisotropy induced by the pre-strain. However, in 
the special case when n is along a principal direction of B and B, the mean 
energy fluxes are along n. 

Also, we note that at the interface C = 0, the mean energy flux (72.)o (say) 
in the substrate is related to the mean energy flux (72.) o (say) in the layer 
through ^ 

ji~'j{n)o = fi-'J{n)o. (6.10) 

Similarly, for the mean energy densities {S)o and {S)o at the interface C = 
, we find 

~-i J (^)o -1 T (^)o rfi m 

^ '{v/vr.y + cos{2J^h) ^ \ + cos{2^hy ^'-^'^ 

Recalling that > -y^, we note in particular that Jl~^J{£)o > fi'^J{S)o. 

We now consider the energy flux velocity defined as the mean energy flux 
vector divided by the mean energy density. For the energy flux velocity g 
wave in the layer, we have 

(g)^ l + coB2e(/, + C) Bn_ 

® {£> (v/v^r + <:oo2J^{h + (n-Bn)' ^' ' 

and for the energy flux velocity g in the substrate, we have 

(n) Bn 
{£) (n • Bn) 

In the layer, the energy flux velocity depends on the depth ( whilst in the 
substrate, the energy flux velocity is the same at all points. This is because 
the wave motion in the substrate consists of a single train of inhomogeneous 
plane waves. On the contrary, the wave motion in the layer may be viewed as 
a superposition of trains of homogeneous plane waves. Because Bn = Bn, 
the energy flux velocities are related through 

[{V/Vr^f + C0S2:^(/1 + C)]^ = [1 + C0S2::^(/1 + C)]^. (6.14) 

Also, because > t^, we note that the energy flux velocity at any point 
of the layer is smaller in magnitude that the energy flux velocity in the 
substrate. Finally we note that 

gn = v, g ■ m = 0, (6.15) 
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in accordance with previous results about finite amplitude inhomogeneous 
plane waves in unbounded deformed Blatz-Ko materials [T7] [TH]. These 
relations are the same as those derived by Hayes in the context of linear 
theories |26j . 

We now define total mean energy densities and total mean energy fluxes 

as 

{s)t= [\m)dc, {n)T= ['mm, 

J-h J~h 

{£)t = / (^(C))dC, (^)t = / mC))dC- (6.16) 
Jo Jo 

The total mean energy density {S)t is the wave energy in the layer {—h < 
( = m ■ X < 0) per unit length (along n) and per unit width (along a) 
of the layer. Similarly, the total mean energy density {S)t is the energy in 
the substrate {0 < ( = m ■ x < oo) per unit length (along n) and per unit 
width (along a) of the substrate. The total mean energy fiux (72.) ^ is the 
energy fiux characterizing the rate at which energy fiows through a normal 
section of the layer {—h<( = m-x<0) per unit width of this section. 
Similarly, the total mean energy fiux {TZ)t is the energy fiux characterizing 
the rate at which energy fiows through a normal section of the substrate 
{—h < ( = m ■ X < 0) per unit width of this section. 
For the wave motion in the layer, we obtain 

{S)r = ipk'a^[lPj^sm{2J^h)+v'h 

{TVjT = ^vPa^ sin(2#^/i) + h\ JlJ~^Bn, (6.17) 
and for the wave motion in the substrate, we obtain 

{£)t = \pk^a\l^ cos^^h), {n)T = lk^a%^-^ cos\^h)ixJ-^Bn. 

(6.18) 

Here we note that the repartition of energy between the layer and the sub- 
strate depends on the depth h of the layer, or more precisely, on the di- 
mensionless parameter kh characterizing the ratio of the layer depth to the 
wavelength. 



7 Interface in a principal plane 

For a given static strain B in the layer and a given unit vector m, there is, 
in general, only one direction n in the interface m ■ x = m ■ x = along 
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which a finite amphtude Love wave as described in Section [S] may propagate. 
Indeed, because n m = n- Bm = is required, n must be along m x Bm. 
However, if m is along a principal axis of B, then n ■ m = n ■ Bm = 
is satisfied automatically for any propagation direction n orthogonal to m 
that is, n can be along any direction in the interface. Here, we consider 
this special case and give a numerical example showing the effects of strain- 
induced anisotropy on the wave characteristics. 

Calling k the unit vectors along the principal axes of B, we write 

n = cos + sin , a = — sinOi + cosOj , m = k, (7-1) 

where the angle 6' G [0, 2tt] is arbitrary. The left Cauchy-Green strain tensor 
in the layer is 

B = Xli®i + Xlj ®j + Agfc ® fc, (7.2) 

where Ai, A2, A3 are the principal stretches in the layer. The left Cauchy- 
Green strain tensor B in the substrate is then uniquely determined as ex- 
plained in Section m First J is determined from equation (14.91) . which here 
reads _ _ _ _ _ 

iG"( J) + fxXfXfJ = lG'{J) + JlX^^X^^Xs, (7.3) 

with J = A1A2A3. Then, using (I4.19p . we obtain 

B = Xli®i + Xlj (g) j + Agfc (g) fc, (7.4) 

where A3 is given by _ _ 

A^ = XfX^'j'. (7.5) 

Hence, the dispersion equation relating the wave speed v and the wave 
number k is (15.111) . or, equivalently, (I5.12p . where 

^J^ = vllc^ = XI cos^ e + Xl sin^ ^, ^J{7^Xl) = vl/{^Xl) = 1. (7.6) 

We note that when both the substrate and the layer are of the Levinson 
and Burgess type (12.71) with Lame parameters A, fi, and A, Jl, respectively, 
the equation (17.31) for the determination of J reduces to the linear equation 

iX + ix + iiXfXf)J-{X + 2fj.) = (X + il + ilXfXf)XiX2X3-(X + 2il). (7.7) 

We now present a numerical example. We take both the layer and the 
substrate to be of the Levinson and Burgess type (I2.4l) - (l2.7p . with X = ^ and 
X = Jl, an assumption often encountered in the geophysics literature (it leads 
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to an infinitesimal Poisson ratio of 1/4, wliicli is common for rocks). Hence, 
G"(J) = 2/x(2J - 3) and G'{J) = 2/1(2 J - 3) and equation ([73D yields 

J = (/I//i)AiA2A3 - 3(/I//i - 1)(2 + K'Xf)-'. (7.8) 

For the ratios po/po and p//i, we take the values, 

Po/po = 1.0, p/p = 0.6. (7.9) 

For the principal stretches in the layer we take 

Ai = 1.45, As = 1.05, A3 = 0.75, (7.10) 

so that J = 1.14, which means a change in volume of 14% . The correspond- 
ing stress tensor in the layer is 

T = p{1.13i ® i + 0.25j ® j - 0.22fc O fc), (7.11) 

so that the deformation can be maintained with the constant normal pressure 
T = —0.22pk (deadload) applied at the upper face of the layer. It then 
follows from (JTS]) and ([73]) that J = 1.18 and 

Ai = 1.45, A2 = 1.05, A3 = 0.77 (7.12) 

As explained in Section [5], the dispersion curves can be deduced from dis- 
persion curves in the linear isotropic case, and are shown in Fig|2l Clearly, 
because Vn varies with the chosen direction n, this figure shows that the num- 
ber of possible modes for a given value of the dispersion parameter kh is not 
necessarily the same for all n. We here focus on the influence of pre-strain 
and choose, for instance, kh = it (wavelength equal to twice the thickness 
layer), a value of this parameter such that two modes of propagation are pos- 
sible in all directions n : a fundamental mode with speed Vi, and a second 
mode with speed V2. In FigOl we plot the polar graphs of Vn, Vn, and of 
the Love wave speeds vi and f 2 of the two modes as a function of the angle 
6 between the propagation direction n and the principal direction i, corre- 
sponding to the greatest stretch Ai = Ai = 1.45. We note that the greatest 
and least values of Vi and V2 correspond to propagation along the directions 
of greatest and least stretch in the interface, indicating an experimental way 
of determining these directions. 

8 Conclusion 

In this paper, we have obtained an exact finite-amplitude Love wave solution 
for a layer and a subtrate consisting of pre-strained compressible neo-Hookean 
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Figure 3: Polar graph of the Love wave speeds vi and V2 corresponding to 
kh = n as a. function of the angle 6 that the propagation direction n makes 
with the direction i. 

materials. The dispersion relation is similar to that of linear isotropic elastic- 
ity, but an explicit dependence on the pre-strain is exhibited. In particular, 
the number of wave modes for different values of the dispersion parameter 
kh {k: wave number, h: thickness of the deformed layer) is influenced by the 
pre-strain. 

It should be emphasized that the existence of this Love wave solution is 
subjected to the condition that the propagation direction n and the normal 
m to the interface are such that nBm = nBm = 0, where B and B are the 
left Cauchy-Green strain tensors characterizing the pre-strain of the subtrate 
and of the layer, respectively. Note that it follows from the continuity of the 
traction at the interface that the conditions nBm = and nBm = are 
equivalent. 

Thus, the interface may be arbitrarily chosen. However, when it is not a 
principal plane of B and B, there is only one propagation direction satisfying 
these conditions. In contrast, when the interface is a principal plane of B 
and B, all propagation directions n in this interface are possible. In this 
case, for a given value of kh, the number of possible modes is not necessarily 
the same for all propagation directions. 

The energy flux and energy density of the solution in the layer and in 
the substrate have been studied in detail. In particular, it has been shown 
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that the mean energy fluxes in the layer and in the substrate are both along 
Bn = Bn, thus along the same direction, parallel to the interface. The fact 
that this direction is not along the propagation direction (except when n is 
a principal direction) is due to the anisotropy induced by the pre-strain. 
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A Uniqueness for the determination of J in 
the substrate 

Here we consider the equation fl4.9p for the detemination of J in the substrate 
when the strain tensor B in the layer is given. It may also be written as 

/(J) = |G"(J) + fiJJ-^m ■ Bm = m ■ Tm, (A.l) 

where T is the stress tensor in the layer, corresponding to the strain tensor 
B. 

First, using the strong ellipticity conditions (12. 6p . we note that 

f'{J) = ^G"{J) + fiJ-^m ■ Bm > 0, (A.2) 



so that f{J) is strictly monotonous increasing for J G [0, oo]. Then, recalling 

lG'{J) < -/i for J < 1, and i( 



G'(l) = -2/i, we have ^G'{J) < ~fi for J < 1, and iG"(J) > -/x for J > 1, 



hence 

lim /(J) = i lim G'{J) < -/i, lim /(J) = oo. (A.3) 

J^O J^O J— >oo 

Owing to the monotonicity of f{J), the limit for J — exists and is either 
finite and negative or — oo. 

If I limj_,o G'( J) = —To < —fi < (with Tq finite), then, clearly, for any 
strain tensor B such that m • Tm > —Tq, equation (lA.ip has exactly one 



solution for J > 0. However, for any strain tensor B such that m ■ Tm < 
—Tq, this equation has no solution for J > 0. ^ 

If limj^QG'{J) = — oo, then, clearly, whatever be the value of m ■ Tm, 
equation (lA.ip has exactly one solution for J > 0. 
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